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Within the Tsallis thermodynamics' framework, and using scaling properties of the entropy, we 
derive a generalization of the Gibbs-Duhem equation. The analysis suggests a transformation of 
variables that allows standard thermodynamics to be recovered. Moreover, we also generalize Ein- 
stein's formula for the probability of a fluctuation to occur by means of the maximum statistical 
entropy method. The use of the proposed transformation of variables also shows that fluctuations 
within Tsallis statistics can be mapped to those of standard statistical mechanics. 



During the last few years there has been a great deal of 
interest in studying nonextensive thermodynamics [Q-^ . 
This results from the assumption of non-additive sta- 
tistical entropies and the maximum statistical entropy 
principle, following the information theory formulation 
of statistical mechanics proposed by Jaynes In- 
deed, besides its relevance in many nonequilibrium prob- 
lems, non-extensivity is of interest for systems of parti- 
cles which show up long-range interactions , as occurs 
in a number of ferroic materials |^ such as ferromag- 
netic, ferroelastic, ferroelectric solids, and astrophysical 
systems H. Within this framework, Tsallis statistical 
entropy l8[ has proven to be the only non-additive gen- 
eralization of Gibbs-Shannon entropy which satisfies the 
following properties: (i) positivity (it takes zero value 
for absolute certainty) increasing monotonously with in- 
creasing uncertainty, and (ii) concavity. However, many 
fundamental features regarding the connection between 
the formulation of statistical mechanics and thermody- 
namics remain unclear. For instance, the identification 
of adequate generalized thermodynamic forces and the 
computation of statistical fluctuations are still contro- 
versial Q. In this letter we clarify such problems and 
provide robust arguments showing the equivalence of the 
present formulations of non-extensive Tsallis thermody- 
namics with the standard extensive equilibrium formula- 
tion. 

Within the Tsallis formalism, the lack of information 
associated with any probability distribution {p{i)} de- 
fined on a set of microstates f2 = {«} is given by 



(1) 



where the parameter q, determining the degree of non- 
extensivity, is positive in order to ensure the concavity 
of <S. The logarithmic function is defined as: lug / = 
(ji-9 _ i)/(i _ g). In the g ^ 1 limit Eq. |l| reduces to 
Gibbs-Shannon's entropy iS = —^^p{i)^T^p{i)- In order 
to simplify the notation, we will indicate the set only 
when necessary. The novelty of the statistical entropy 
(|l|) is that it does not satisfy additivity. Instead, for two 
systems A and B described by independent probability 
distributions: 



S{A + B)= S{A) + S{B) + (1 - q)S{A)S{B). (2) 

For a given physical system the equilibrium probabil- 
ity distribution {p*{i)} corresponds to the distribution 
that maximizes S under the normalization condition 
(Eier2P(*) ~ 1) ^^'^ relevant constraints imposed 
by the available statistical information on the system. 
The thermodynamic equilibrium entropy is then iden- 
tified with kS*, where S* = S{{p*{i)}) and k is the 
Boltzmann constant. In the case of an isolated (micro- 
canonical) system no statistical information is available 
and maximization of (|l|) leads to the equiprobability dis- 
tribution on fl. Difficulties arise when trying to study 
non-isolated systems for which some constants of motion 
Xa{i) in the isolated system, such as the energy, the vol- 
ume, the magnetization, the number of particles etc., 
become fixed only on average ((Xa)). The index a ex- 
tends over the number of such controlled observables of 
the non-isolated system. It is worth noting that a crucial 
property of such observables, required by the foundations 
of thermodynamics |lO[ , is that they are additive, i.e. for 
two independent systems A and B: 



{X^{A + B)) = {X^{A)) + {X„{B)) 



(3) 



Within the Jaynes scheme, imposing knowledge of the 
statistical information (i.e. the controlled average values 
{Xa)) on the maximization procedure, leads to the stan- 
dard, canonical, isobaric, isofield, grand-canonical, etc... 
ensembles. Two different strategies can be adopted when 
using TsalHs statistical entropy ^: 

(i) Unbiased averaging: when the average quantities 
of the microscopic properties are defined in the standard 
way {Xa) = X^iGO Legendre structure of 
standard thermodynamics is not recovered i.e. in 
equilibrium, the partition function cannot be related to 
a Legendre transformation of the entropy. For this reason 
this scheme has not been studied very much. 

(ii) Biased averaging: A Legendre structure is 
recovered if averages are defined as: {Xa)q = 
'^i^QPq(,'i)Xa{i), where the averages are performed 
using the so-called q-scort probabilities: Pq{i) = 
[p(i)]^/ Xjgj2[p(i)]*. After maximization when this 
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choice is used, one recovers a fundamental thermody- 
namic identity which reads: 



dS* = ^yad{Xa)q, 



(4) 



where ya correspond to the Lagrange parameters nec- 
essary for keeping the average constraints {Xa)q when 
performing the maximization of the statistical entropy. 
The physical meaning of the parameters ya is, at present, 
still controversial , and they are not necessarily those 
that control mutual equilibrium between thermodynamic 
systems (physical temperature, pressure, magnetic field, 
...). Besides, the equilibrium statistical entropy can be 
expressed, in any statistical ensemble, as S* = lnq Zq, 
where Zq plays the role of the partition function and can 
be written as: 



(5) 



with 



= [X ~ (X),)/E,enb«]'- . The q- 
exponential function is the inverse of In^ and is defined 

[1 + (1 — q)f]^^ ■ This definition only makes 



as 



sense for 1 + (1 — q)f > and, in practice, this means 
that within this biased scheme only the range < g < 1 
has a physical meaning p^ , p^ . 

It is worth noting that a third "averaging" procedure, 
consisting of using the above proposed biased averag- 
ing, but without normalization, has been shown |^ to 
be equivalent (after suitable redefinition of the Lagrange 
multipliers) to this normalized biased averaging strategy. 

At this point it is convenient to analyze the scaling 
properties of the non-extensive statistical entropy S sat- 
isfying the property (|^). By recursively applying equa- 
tion (0) to A identical systems one obtains: 



5(AA) = ) 5(A) + (^) (l-q)5(A)2 + ... 
■■■ + {^)il-q)'-'SiA)\ 
which can be easily expressed in the compact form: 



S{XA) 



1 
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{[l + {l-q)S{Ar~l} 



(6) 



(7) 



Of course, for g ^ 1 the usual scaling for extensive sys- 
tems is recovered. 

Let us now concentrate on the analysis of the biased 
scheme, and consider the statistical entropy iS* of the 
equilibrium states as a function of the set of extensive 
quantities {Xa)q- Under the assumption that expression 
(^ [and therefore (|^)] can be applied to the equilibrium 
distributions and thus to 5*, we can substitute S*{A) 
by S*{{{Xa)q})- Differentiating the expression (|^) with 
respect to A and taking A = 1 we obtain: 

Y,{Xa)qyo. = [1 + (1 - ci)S*] In [1 + (1 - q)S*] . 

a. ^ 

(8) 



Equation (^ represents a generalization of Euler's theo- 
rem for the Tsallis non-extensive entropy. By differenti- 
ating (||) and using the thermodynamic identity (||), we 
deduce that: 



J2iXc.)qdya = In [1 + (1 - q)S*] dS* 



(9) 



which corresponds to the Gibbs-Duhem equation for non- 
extensive systems. It indicates that the Lagrange param- 
eters ya are not intensive variables except for q = 1. 

The equation obtained suggests a transformation of 
variables that keeps the structure of the thermodynamic 
identity and allows to recover the standard Gibbs-Duhem 
equation. Indeed, defining: 



S* 



dS* 



l + {l-q)S* 



In [1 + (1 - q)S* 
1-q 



(10) 



and 



Va 



"'^ l + (l-g)5*' 
expressions (^ and (|^) can be rewritten in the form 

dS* = y^^yad(Xa)q, 



(11) 



'^{Xa)qdya = 0. 



(12) 



(13) 



Therefore, iS* is an extensive statistical entropy (satis- 
fying 5*(AA) = XS*{A)) and ija are the conjugated in- 
tensive variables of {Xa)q- It should also be noted that 
the condition q < 1 (necessary for the transformation to 
be well defined) ensures that S* and S* have the same 
concavity (thus the same maxima and minima). 

The above transformation of variables (^0|) and (ll|) 
has partially been suggested by different authors dur- 
ing recent years. Tsallis has already shown the extensive 
properties of 5* , which in fact had been proposed sev- 
eral years before by Renyi |14|. The intensive character 
of the variables ija was recently pointed out by Abe et al. 
|]l2| , who showed that such variables control the mutual 
equilibrium between thermodynamic systems. Finally, 
the fact that standard thermodynamics is fully recovered 
after the transformation of variables proposed above has 
also been implicitly suggested by Toral ||l^ within the 
framework of the microcanonical ensemble. 

Another aspect which is still under debate refers to 
the calculation of statistical fluctuations within Tsallis 
statistics. The computation of the variances of extensive 
quantities {X"^) has not been clearly solved ||l^ not only 
within the unbiased averaging scheme, but also using bi- 
ased averaging. Here we provide a method of computa- 
tion of fiuctuations which is based on a generalization of 
Einstein's fluctuation formula and is independent of the 
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averaging scheme. Moreover, for the case of biased aver- 
aging, the use of the transformation of variables proposed 
above, allows standard thermodynamic fluctuations to be 
recovered. The method has been recently proposed for 
the study of fluctuations in standard extensive statistics 

In order to compute fluctuations within a very gen- 
eral framework we must study how the distribution p*{i) 
changes under virtual displacements. Such displacements 
can be understood as a consequence of an internal con- 
straint which causes the deviation of the system from 
the equilibrium state, or as a consequence of a contact 
with an external bath which allows changes in parame- 
ters which would be constant under total isolation. Let us 
consider that such a displacement can be characterized 
by a multidimensional parameter k which takes values 
on a continuous set TZ and vanishes at equilibrium. We 
define the set of functions {p*{i\K)} as the displaced (con- 
ditional) probability distributions for a given k obtained 
after maximization of (|l|). Note that p*{i\0) = p*{i)- 

The key point for characterizing fluctuations is to con- 
sider the joint probability distribution p{i, k) defined on 
the enlarged set x TZ. This means that we are assum- 
ing that K is a fluctuating variable and that the system 
is described by simultaneously specifying i and k. The 
problem now is to determine p{i, k). Using the same phi- 
losophy used for the determination of equilibrium dis- 
tributions p*{i\K), we propose to generalize the maxi- 
mum statistical entropy principle to the determination 
of p*(i,K) on 17 X 7^. This means maximizing: 



SnxTZ {{p{i,K)}) = - 



TZ 



[p{i,K)]'^lnqp{i,K)dK, (14) 



where we have assumed that the measure of the lack of 
information on 17 x 7?, is also given by Tsallis entropy . 

Taking into account the mathematical definition of 
the conditional probability p(i, k) = p{K)p*{i\K) and us- 
ing the standard property In, fg = lug / -t- lug 5 + (1 — 
q) lug /lug o and the normalization condition for p*{i\k), 
equation (y^ can be written as: 



[p{^WS*^i^)+STz{{p{K)}), 



(15) 
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where S^{k) — Sfi ({p*(«|k)}) is nothing more than the 
thermodynamic entropy of the virtually displaced sys- 
tem. Therefore, the problem of determining the distri- 
bution of fluctuations reduces to the maximization of the 
functional ( [T5| ) with respect to p{k) under the normal- 
ization condition J^p{K)dK = 1. This yields: 



In 



(16) 



Since the function converges to the standard exponen- 
tial function in the limit g — > 1, the above equation just 



reduces to the well-known Einstein formula in the exten- 
sive limit. This probability distribution for fluctuations 
P*{k) is, therefore, the most probable (or the least bi- 
ased) distribution of the fluctuating parameter k. Also 
note that formula ( p^ can be deduced in a non-rigorous 
way (as done in many statistical mechanics text-books in 
the case of extensive statistics) by inversion of the gener- 
alized Boltzmann formula for the microcanonical ensem- 
ble: S* 

Taking into account the fact that S*{k) is maximum 
at K = we can expand it in powers of k, as: 



luq W , where W is the number of states in 



«2 + 



(17) 



with {d'^S*/dn'^) < 0. Substituting in (|T^), we find that 
small K fluctuations are Gaussian distributed according 
to: 



P*{k) ~ exp { -- 



1 

2 r 



d-'S' 



{l-q)S*{0) 



(18) 



This formula, as well as (16), is independent of the differ- 
ent averaging procedures explained in the introduction. 
It shows that the probability distribution for fluctuations 
of measurable quantities depends on q as well as on the 
entropy of the reference equilibrium state around which 
fluctuations take place. 

In order to compute the variance (k^) the two aver- 
aging schemes discussed above can be considered. Nev- 
ertheless, if K is a macroscopic quantity, the unbiased 
scheme seems more appropriate. The two choices yield 
very similar results in the limit of small fluctuations. If 
we use unbiased averaging, the variance is given by: 



[{l + {l-q)S*m 



(19) 



whereas in the case when we use the biased averaging 
procedure, one gets the same result, but it is divided 
by q (the q-scort probability of a Gaussian distribution 
is also Gaussian). It is worth noting that within the 
two interpretations, as g < 1 fluctuations are larger than 
those corresponding to the extensive case 9=1. 

It is also interesting to note that the transformation of 
variables ( p^ ) and (|Tl|) suggested in the previous section 
allows standard thermodynamics to be recovered within 
the biased scheme, and can be used to rewrite equation 
( |T6| ) in terms of S* . Using the fact that S* and S* have 
the same concavity, a straightforward calculation yields: 



P*{k) 



(20) 



which is the standard Einstein formula, and the variance 
transforms into: 
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(21) 



Therefore the transformation of variables aUows not only 
standard thermodynamics to be recovered, but also the 
same statistical fluctuations. For instance, a simple cal- 
culation for the energy fluctuations AE from this expres- 
sion gives: 



{AE-^) = - ( ^ 

^ dp 



(22) 



where = f3/[l -I- (1 — q)S*] is the intensive physical 
inverse temperature, and (3 is the Lagrange parameter 
associated with {E)q. 

In summary, we propose that Tsallis statistics, within 
the biased averaging scheme, can be mapped into stan- 
dard thermodynamics. This is done by a suitable trans- 
formation of variables, which appears after the deriva- 
tion of the generalized Gibbs-Duhem equation resulting 
from the scaling properties of Tsallis entropy. The inten- 
sive generalized forces in such framework arise naturally 
from the transformation of variables. In addition, the 
standard equilibrium fluctuation-dissipation theorem is 
also recovered. This indicates that, besides the thermo- 
dynamic equivalence, statistical fluctuations also behave 
according to standard statistical mechanics. Moreover, 
this ensures the equivalence of the equilibrium response 
functions. 
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